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Abstract 

We study the finite-temperature phase transition of the generalized Gross-Neveu model 
with continous chiral symmetry in 2 < d < 4 euclidean dimensions. The critical exponents 
are computed to the leading order in the 1/N expansion at both zero and finite temperatures. 
A dimensionally reduced theory is obtained after the introduction of thermal counterterms 
necessary to cancel thermal divergences that arise in the limit of high temperature. Although 
at zero temperature we have an infinitely and continously degenerate vacuum state, we show 
that at finite temperature this degeneracy is discrete and, depending on the values of the 
bare parameters, we may have either total or partial restoration of symmetry. Finally we 
determine the universality class of the reduced theory by a simple analysis of the infrared 
structure of thermodynamic quantities computed using the reduced action as starting point. 
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1 Introduction 



One of the differences between field theories and mechanical systems is that field theories have 
an infinite number of degrees of freedom, which makes possible spontaneous symmetry breaking. 
More generally, we know, by now, that also dynamical symmetry breaking plays an important 
role in modern physics. In condensed matter physics, for example, such a mechanism is used to 
describe superconductivity, by the condensate of Cooper pairs, while for particle physics it is the 
main source of hadron masses and governs the low energy hadron dynamics. 

Connected with dynamical symmetry breaking comes the problem of restoration of broken 
symmetries for some sufficiently high temperature. As it is well known, a phase transition occurs 
when there is a singularity in the free energy or one of its derivatives. The fact that almost 
all phenomena studied in theories near the transition point exhibit scaling, that is, a power- 
law behavior between two mesurable quantities, leads naturally to the classification of different 
conformal field theories in universality classes, which are itself determined by a set of numbers 
usually called critical indices. 

A large number of conformal field theories is known in two dimensions and none is firmly 
established in four. The dimensionality (d = 3) is between these extremes and contains two 
well established families of conformal field theories. The first family contains the usual Ising, 
XY and Heisenberg critical a models and is characterized by a particular set of exponents called 
scalar critical exponents. The second family contains various critical four fermion models Jl| 
and the set of exponents that characterizes it is called chiral critical exponents. The critical 
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properties of the great majority of phase transitions in three euclidean dimensions (magnetic 
systems, superconductors etc.) are quite accurately described by the first family of conformal field 
theories while the finite-temperature phase transitions in certain (3 + 1) dimensional quantum field 
theories are argued to belong to these same universality classes 0. 

When studying the finite-temperature chiral restoration in QCD one is usually gided by the 
concepts of dimensional reduction and universality. From the dimensional reduction point of view, 
a hot field theory can be regarded as a static field theory at zero temperature in a lower dimension 
0. For example, four dimensional QCD with Nf light quarks near the transition can be described 
by the three dimensional linear a model with the same global symmetry M. On the other hand, 
as the classification in universality classes is done by the computation of critical indices and those 
critical indices are infrared (IR) sensitive, we must understand the role of each degree of freedom 
in this context. This can be easily done with the use of effective field theory methods in which the 
integration over different energy scales gives a theory for the degrees of freedom which are the real 
responsible for infrared divergencies ||. Dimensional reduction is based in the zero temperature 
Appelquist-Carazone decoupling theorem ||, where a low energy theory is constructed by the 
integration over heavy fields in the functional integral. 

According to standard dimensional reduction arguments, the fermions themselves, even if 
they are massless at zero temperature, do not influence the nature of phase transition at finite- 
temperature. It is rather their bosonic composites, the Goldstone bosons, which are of importance. 
This follows directly from the universality of second order phase transitions [[!]] , in which the com- 



2 



monly held assumption is that all the possible universality classes (or equivalently, conformal field 
theories) are variations of the a model and one need only match the correct symmetry-breaking 
pattern. As a consequence, the chiral transition of four dimensional QCD, with Nf = 2 flavors, 
should lie in the same universality class as a three dimensional 0(4) magnet. Similarly, other 
models, e.g. four-fermion theories in d dimensions such as the Gross-Neveu model (GN model) 
H with discrete symmetries or the Nambu-Jona-Lasinio model (NJL model) || with continous 
chiral symmetry, are expected to be in the same universality class of a (d — 1) dimensional Ising 
or Heisenberg magnet, respectively. 

Recently, however, it was pointed out that there exist different (d = 3) conformal field theories 
with the same symmetry-breaking pattern Q. They are exactly the four-fermion interaction 
models of the Nambu-Jona-Lasinio type. Physically, this corresponds to the fact that on the 
chirally symmetric side of the phase transition there are N massless fermions whose effect is 
felt even in the IR fixed point, just like the Goldstone bosons. The presence of more than one 
universality class in (d = 3) makes the procedure of dimensionally reducing a quantum field theory 
ambigous and it is now uncertain to which conformal field theory the (d + 1) finite-temperature 
quantum field theory will reduce. The argument if favor of the bosonic universality class goes as 
follows. At finite-temperature, the (d = 4) fermion reduces to a collection of (d = 3) massive 
fermions and there is no zero mode for which the Matsubara frequency vanishes. Nevertheless, 
even if a single massive field does not influence the phase transition, the cumulative effects of an 
infinite number of such fields may have an appreciable impact. In order to see whether or not this 
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happens, all harmonics should be summed and their cumulative effects studied. 

It is the purpose of this paper to discuss the assumptions underlying this analysis and to 
determine to which universality class the simplest generalization of the Gross-Neveu model with 
continous chiral symmetry, belongs. Besides being an interesting theoretical model, it is also 
believed that, when properly extended to incorporate continous chiral symmetry, four-fermion 
models are more realistic as effective theories of QCD than the linear a model, especially at scales 
where quark structure is important. 

The paper is organized as follows. In section II, we present the model. In section III we com- 
pute the critical exponents at zero temperature while, in section IV, we obtain, after dimensional 
reduction, the new critical exponents that determines the universality class of the reduced theory. 
In section V, we analyse the vacuum structure at finite-temperature. Conclusions are given in sec- 
tion VII. In the appendix, we compute the thermal renormalization group functions that controlls 
thedependence of the thermal counterterms on the temperature. In this paper we use ^ = c = 1. 

2 The chiral-continous GN model 

We are interested in studying the behavior of a multiplet of N fermions coupled with a pair of 
composite self-interacting pseudo-scalar and scalar fields, cr(x) and n(x) respectively, in such a 
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way that the euclidean functional action reads 

Se = J d d * {-</; + g B (a + i ls *)) V» + ^(<V) 2 + i(9 M 7r) 2 - im| (a 2 + tt 2 ) + ^ (a 2 + vr 2 ) 2 } 

(1) 

This model has continous chiral symmetry 

and <7 — > — cr, 7r — > 7r with a being a constant and we use for the 7 S the following representation 

In order to have an uniquely defined vacuum (see fig. 1), we must choose a prefered direction 
in the space of components, say (a) = v and (ir) = 0, and shift our fields as 



(2) 



I I ^ 



-1 



cr(x) — > v + <j(x) 
ir(x) — > n(x) 

where v = In this new picture the euclidean action becomes 



(3) 



Se = /^{-im|t; 2 + ^; 4 -^ + ^ (4) 
- \mlo 2 - \m 2 y - X B va 3 - X B vaTr 2 + ^ (a 2 + tt 2 ) 2 1 , (5) 

with m CT = m 2 B — 3\ B v 2 and = m 2 B — A^f 2 . If we use the definition of v we can see that 
m n = 0, as it should be since the n field is a Goldstone boson. Although new interaction terms 
have arised, it is easy to prove that this renormalizable model needs only usual wave function, 
mass and coupling constant renormalization to render the theory finite. 



Figure 1: Mexican hat shape for the zero temperature effective potential in the tree level. The vacuum is 
infinitely degenerate and one of the degrees of freedom will become a Goldstone boson after the definition of 
an unique ground-state. 

The finite-temperature version of this model is obtained, in the Matsubara formalism, by the 
compactification of the imaginary time dimension together with the imposition of (anti)-periodic 
boundary conditions to (fermionic)-bosonic fields as 

^(x,r) = -= J2 &(x)e** T , (6) 

VP n=-oo 

with ui n — jf(n + |) and 

-| oo 

°^ T ) = ^n E °n{*)e^\ (7) 

VP n=-oo 



with ui n = 

Since global aspects of the euclidean manifold do not affect local properties of a quantum field 
theory we will, again, need only those previously mentioned renormalization constants to render 
the theory finite. 



3 Critical exponents at zero temperature 

The chiral-continous Gross-Neveu model, in contrast to the discrete Gross-Neveu model, is renor- 
malizable in four dimensions and it can be described, beyound the tree level, by all the techniques 
developed for the ($ 2 ) 2 theory: renormalization group equations and large N expansion, for ex- 
ample. 

To investigate the vacuum structure of the theory, we must consider the leading order density 
V e ff evaluated for constant fields (a) = v and (n) = ]TD[. This gives the gap equation (a cut-off 
A is implied) 

^ = ^_ mSj + w = 29l | A |^_L_. (8) 

We define the critical coupling g c as 

9 , r A d d q 1 



so that the gap equation takes the form 



ILh-^\- f A ddq gW do) 

2 { mBt t + l)~J (2iryq%q*+g%v^ 



where t = f — lj is the reduced coupling. 

This form for the gap equation is particularly well suited for extracting critical indices since the 



problem reduces to counting the infrared divergences on the right hand side [II]. As previously 
said, the critical indices are the responsible for the power-law behavior between two mesurable 
quantities as we approach the critical coupling. In this sense we will compute the critical exponents 
defined by (vpipj ~ t 13 , ~ v 1 ^ 5 , (S)|^ ~ t" etc, and show that the introduction of 

a self-interaction for the bosonic fields does not change the values of the exponents obtained for 
the NJL model. 

First, we can see that the critical indices for the NJL model are recovered by simply putting 
X B = in the gap equation (10). Indeed, in this case the infrared behavior of the integral with v 



goes as 

f A g B v2 „2„,d-2 4 d — 4 6 d 6 

/ in \d 2t 2 ; — 2~^\~9bv - g B v + g B v (ii) 

and the critical indices are simply (3 = u = \/{d — 2) and 5 = d — 1, since, to the order we are 
concerned, the self-energy and the two point function coincide. 

Now, let us show that even with A ^ we still have the same critical indices. For d < 4 this 
is rather obvious since the singularity of the integral in (|10D allow us to neglect the v 2 term in the 



left hand side. However, if we are in four dimensions, we can use the integral behavior (|TT|) to 
obtain 

2 9 2 B v 2 , (12) 
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and compute the f3 exponent from the behavior of 



m 2 D t 



v2 ^ (13) 
TTi+ 2 9% 



with t. This gives the mean field value /3 = |, a characteristic of conformal field theories in 
dimensions greater or equal to 4. 

We conclude that the introduction of self-interaction for the composite bosons, has not changed 
the structure of the IR singularities. Also, we can see that, in (d = 4) all critical exponents are 
mean field and there are no non-gaussian fixed points while, in three dimensions the value (3 = 1 
tell us this is a chiral conformal field theory of the Nambu-Jona-Lasinio type. 

4 Dimensional reduction and universality 

We will now consider the problem of computing critical indices when a non zero temperature is 
introduced. We could do this in a similar way to what we have done in the zero temperature case 
but we will, instead, compute the fermionic determinant using a modified minimal subtraction 
scheme at zero momentum and temperature 1/(3 because this is a more suitable way of obtaining 
the functional form of the dimensionally reduced theory. 

As we have pointed out in the introduction, the reduced theory is an effective theory for the 
zero modes of the original fields || (and see also 10). It can be explicitly obtained from the 
initial partition funcion by integrating out the uo n ^ Matsubara frequencies of the bosonic fields 
and, as we do not have a zero mode in the case of fermions, the complete integration over the ip 
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and ip fields. The resulting QFT will be described by an euclidean partition function of the type 

Z = J Va V7i exp{- J cf^x (£ e //K tt ) + 5C) }, (14) 

where 5C includes all other local terms not present in the original lagrangian density but which 
are consistent with the symmetries and, in the lagrangian density £ e //(o" , 7r ), the parameters will 
be, in general, functions of the cut-off, the temperature and the bare parameters. 

Since we are interested in the computation of the finite-temperature effective action at the 
leading order in the 1/N expansion it is convenient to define 

9 = 9b^N 

< m 2 = m 2 B N , (15) 
A = X B N. 

and 

rhi = g B (cr + m Q ) 

(16) 

m 2 = g B (oo - ino) , 
being the eigenvalues of M = (?_b(I(To + ils^o), where I is the unit matrix. 

At finite-temperature, the integration over the imaginary time becomes a sum over Matsubara 
frequencies and if we remember that the fields ip and tp are anti-periodic in the time component, 
we obtain, after integrating out the fermionic degrees of freedom, the expression for the fermionic 
determinant 

o oo . id— It, 2 

trln( 9 + M) = i Y.J JLJL gl„ + ^ + *J), (17) 
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where, as usual, u; n = ^(n + |) for fermions, or yet 



O CO „ 

trln(cj) + hi) = - J2 J 



2' 



(27T) 



d-1 



2 In f k 2 + ujI 



In 1 + 



k 2 + u; 2 



Zn 1 



mo 



k 2 +u; 2 



The first logarithim in the above expression is naturally absorved as an overall normalization 
factor. The other two can be written, after integration over k and summation over n, as a series 

:-y~g 2s 



trln($ + M) =2J2 



s=l 



(47r)^ r( s ) 







d-2s 



(19) 



where £ is the analytic extension of the modified Epstein- Hurwitz zeta function [13]. 

In the limit of high temperature (/3 — > 0) the only contributions to the effective action comming 
from expression ([II]) are those given by s = 1 and s = 2. Those contributions are divergent so that 



thermal couterterms will be needed to make the reduced theory finite in this limit |jl2| . In this 
context, we join the above divergences with those comming from the integration over non-static 



bosonic modes [ T4j and introduce the following counterterms (for d — 4 — e) 



5m2 = ohM A ~ XI ~^ s 



Am 2 I (I , 2/3u 

i^u +i --'- ta - 



and 



10 A 2 g A 

T _|_ J 

6 64tt 2 4 16vr 2 



1 -3 
-3 1 



7T 



1,7,, 2/V 
h - + In- 



(20) 



7T 



(21) 



where, as usual, 7 is the Euler number, /1 is a mass parameter introduced by dimensional regular- 
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ization to give Green functions their proper dimensions and we have defined, for simplicity 



1 1 



(22) 



Since as (3 — > we are in the disordered phase (v — 0), we define 





(23) 



and 



(T 



so that the renormalized action for the (d = 4) becomes 



(24) 



S ef f(a ,7r ) 



J> , 2 8A 3 C(3)^ a 







W + 



9 2 9 7T 4 



(25) 



(26) 



where, as usual, 

m 2 ((3) = -m 2 I + Sm 2 
\(/3) = AI 4 + <JA. 

Now we are ready to start the computation of the critical indices for the reduced theory. In 
this sense we will, again, set (n) = and expand our field a around a constant configuration v to 
obtain the thermal effective potential 



(27) 



4/3 9 2 9 tt 4 

The first critical exponent one usualy computes is the (3 exponent which gives the power- 
law behavior of the order parameter with the temperature. The critical temperature f3 c is itself 
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determined by the requirement of a vanishing thermal mass when (3 — > (3 C . In this sense the 
condition m 10 (/3 c ) = gives 

In order to have a critical temperature we will, for the moment, consider the case 4A > g 2 . Now, 
we can rewrite m\ Q {(3) as 

<W = ^-f)^-^j. (29) 
As (3 — > f3 c we see that, by defining the reduced temperature 6 as 

l3 = T+Wc' (30) 

we obtain 

m\M = Y A {^~ 9 2 ) J c + 0(9 2 ) ~ a + b6 : (31) 
that is, m\ Q {pi) is linear in 9, or equivalently, we have for the first critical exponent the value (3 — \. 

The second critical exponent, the 5 critial exponent, is obtained from the inhomogenious gap 
equation, when we consider an external source J for the a field. In this case, the gap equation 
becomes 

m ;„03). + ^». 3 + f^ 5 = J, (32) 

and, as we approach criticality (m w ((3 c ) — > 0), the gap equation gives 5 = 3, since v 3 dominates 
the v 5 term. Finally, as we are in the leading order in the 1/N expansion, once again the two point 
function and the self-energy are the same. This tells that the v critical exponent is also equal to 
(3 assuming the value |. 
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After all, to which universality class does the reduced theory belong? As we have seen, all 
critical exponents are mean field in the leading order of the 1/N expansion. Actually, one would 



get the same critical exponents in this approximation in d = 3 [[L5[| , although they would obviously 
not satisfy hyperscaling relations in this dimension. Notice also that they are almost as far from 
the scalar critical exponents (3 ~ 0.3 in three dimensions J7j, as from the chiral critical exponents 



computed in the previous section (3 = 1 [16 



We would be tempted to conclude that dimensional reduction really spoiled the critical behavior 



of thermodinamic quantities ||17| , reforcing the thesis that this procedure is ambigous. However, 
this sentence must be used with care. The critical exponents we have just computed are related to 
an effective field theory obtained from the integration over all fermions and all nonstatic bosons. 
Since the critical indices come from the IR behavior of mesurable quantities near the transition 
point and the integration we have done gives IR finite results, it is reasonable to expect that the 
above degrees of freedom will not improve mean field values for the exponents. In order to recover 
non trivial values for the critical indices we must consider the effects of the zero modes <To and 
ttq in the computation of thermodynamic quantities, because they are the only responsible for IR 
divergences in the finite-temperature version of this model. 



Following the same steps of Kogut in [|Ty| we will consider the critical behavior of the suscep- 
tibility x f° r a a model with action given by eq. (J^|). Defining the critical curvature \j? c as the 
point where the susceptibility diverges 



> XiM f A d d ~ l q 1 



^~ (3 J (2tt)<*-V ( 33 ) 
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the expression for the inverse susceptibility can be recast into 

_ x / , A 10 ( I) /'A d d -\ 1 , _ , , 



Again, the extraction of the critical index 7 reduces to counting the powers of the infrared singu- 
larities in the left hand side of the previous equation. In three dimensions (d — 1 = 3) the second 
term in eq. ([34] ) dominates the scaling region, giving the zero temperature susceptibility exponent 



7 = 2/(d — 3). It is easy to obtain the other critical exponents and they show the same type of 
behavior as 7. The reduced theory has scalar critical indices and due to the dimensionality of the 
order parameter we conclude that it belongs to the universality class of the Heisenberg magnet, 
as reasonably expected. 



5 Vacuum structure at finite-temperature 

In the previous section, we defined the vacuum state in the o"o direction, breaking spontaneously 
a symmetry, and computed the critical indices from the effective potential (p7|), obtained after 
dimensional reduction. Now comes the question. Is this really a choice? As it will become clear 
soon, this is definitely not a choice. 

If we still had the shape of a "mexican hat" for the effective potential (see fig. 1), we would 
simply have to choose a prefered direction to be the minimum configuration defining the vacuum 
of the theory. However, differently from the zero temperature situation, this is no longer the case. 
As we can see from expression (f20|), the integration over fermi fields give different thermal contri- 
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butions for the masses of o"o and ir and if we define the critical temperature as the temperature 
in which we have a vanishing mass, we will conclude that, since we have two different masses, we 
will have two critical temperatures f3~ l and such that 




Figure 2: Vacuum structure for j3 > (3 n . We see that the ground state is twice degenerate in the ao direction 
while a false vacuum (or metastable vacuum) can be identified in the ttq direction. After the identification of 
the true vacuum state we will no longer have a Goldstone boson. 



1 (4A - g 2 ) 



24 
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(35) 



and 

1 U\ + g 2 ) 

Now, let us show that associated to each critical temperature there is a second order phase 
transition and a symmetry restoring phase. Indeed, if we go from zero temperature {(3 = oo) to 
a temperature close to f3~ l the effective potential takes the form of (fig. 2). There is a discrete 
degeneracy of the vacuum showing that in the previous section we made a correct choice of the 
ground-state. 

If we continue increasing the temperature in such a way that f3 a < f3 < (3 n we note (see fig. 3) 
that the discrete symmetry in the 7To direction is restored but we continue with a broken phase in 
the o"o direction. Complete restoration of symmetry (fig. 4) will only be achieved when (5 < (3 a , 
and this is exactly the picture after dimensional reduction since we have already taken the limit 
0. 

All the above discussion is valid only if 4 A > g 2 . What happens when 4A < g 2 ? In this case 
we will have only the f3 n critical temperature. This is to say, the discrete symmetry on the no 
direction will be restored while, in the o"o direction, the minimum will become deeper and deeper 
as we increase the temperature. We will not have a complete restoration of symmetry and the 
effective potential will have the shape of (fig. 3). 
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Figure 3: Partial restoration of symmetry for f3 a < (3 < (3 W . The ttq degree of freedom is now in the disordered 
phase while we still have discrete breakdown of symmetry in the <7rj direction. 

6 Conclusions 

Dimensional reduction of the chiral-continous Gross-Neveu model was performed giving an effective 
theory that belongs to the same universality class of the Heisenberg magnet. In spite of the 
integration over fermions and nonstatic bosons modes have given mean field results for the critical 
indices, we were able to explicitely determine the universality class of the reduced theory by 
simply taking into account the effects of the only relevant degrees of freedom in the computation 
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Figure 4: Complete disorder for (5 < (3 a . We see that, due to the difference in the thermal masses of the 
static fields <7n and ttq this surface cannot be obtained by the revolution of a generating curve. 

of thermodynamic quantities; the static Matsubara modes. 

The presence of a coupling 275 in the original action is the responsible for the difference on 
the thermal contributions for the masses of the pseudo-scalar and scalar fields a and 7r . This 
difference has important implications. First, we saw that if 4A > g 2 we have restoration of 
symmetry for some critical temperature, first in the 7r direction and after in the a direction, 
while, if 4A < g 2 we have only partial restoration of symmetry. The introduction of temperature 
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in a theory with a continous chiral symmetry has changed the role played by the ttq field from the 
zero temperature case, where it was a Goldstone boson, to the finite-temperature case, where this 
is no longer true. 
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Appendix A - Thermal RG flow 

We compute, for completeness, the thermal renormalization group functions which controlls the 



dependence of the counterterms on the temperature. As it is well known |12|], the transition from 
a four- dimensional theory to an effective three-dimensional one is renormalization point depen- 
dent. This is not surprising since the Appelquist-Carazone decoupling theorem || only holds if a 
particular class of renormalization prescriptions is adopted, being optimal in BPHZ subtractions 
at zero momenta and temperature (3. Neverthless, we can write down the thermal renormalization 
group equation which controls the dependence with the temperature on the correlation functions 
and compute the running quantities. 

The independence of bare correlation functions on the choice of the renormalization point, is 
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expressed by the conditions 



*«f rw = o, 

T J_ T (N) = o 

dT 



(37) 
(38) 



and, for simplicity, we will work, in this section, with the temperature T instead of (3 = 1/T. 

Since we are only interested on the computation of thermal renormalization group functions, 
we will consider solely 

/ f) f) f) \ 

(39) 



where 



(3t = T 



d\(T) 
dT ' 



-2 T dm 2 (T) 



and 



r (JV) 



^eff 



dog 



o-o=0 



Using for the effective potential the expression (|27|) we easily get 



and 



7t 



-f4A 
12 V 



£ 2 



m 2 47r 2 



(40) 
(41) 

(42) 



(43) 



(44) 
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Equation fl39| ) can be solved by the method of characteristics. One simply introduces a dilata- 
tion parameter t and look for functions X(t) (for fixed g) which satisfies 



with A(l) = A. If we now solve the above equation we get 



To investigate the large T limit we will have to study the behavior of the effective coupling 
\{t) as t goes to zero. Since (3(\) is positive \(t) will decrease. Moreover, since the slope of the f3 
function is also positive, the gaussian IR fixed point is attractive. 

We can solve a similar equation to g and defining (A*, g*) as the zero of the f3 function 



we see that, although the anomalous dimension for the <To field may change when we pass from 
the situation of total to the situation of partial restoration of symmetry, the i] exponent defined 
as 7] = 7](X*,g*) remains unchanged (rj = 0). 

One may ask whether renormalization group improves the critical indices after dimensional 
reduction or not. The answer is no and this should not be surprising since the form of the 
ultraviolet and thermal divergences that we had to renormalize are of the same type as for a 
theory in four dimensions, where everything is mean field. 




(46) 



P(\*,g*) = 0, 



(47) 
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